Verifying entanglement in the Hong-Ou-Mandel dip 
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The Hong-Ou-Mandel interference dip is caused by an entangled state, a delocalized bi-photon 
state. We propose a method of detecting this entanglement by utilizing inverse Hong-Ou-Mandel 
interference, while taking into account vacuum and multi-photon contaminations, phase noise, and 
other imperfections. The method uses just linear optics and photodetectors, and for single-mode 
photodetectors we find a lower bound on the amount of entanglement. 



I. INTRODUCTION 

Quantum interference effects that arise when single 
photons impinge on a beam splitter are crucial to linear- 
optics quantum computing schemes [hri^ , with the other 
indispensable nonlinear ingredient provided by photon- 
counting measurements. One such linear-optics quan- 
tum interference effect was observed for the first time in 
1987, by Hong, Ou, and Mandel, and it still carries their 
name j4j. In the Hong-Ou-Mandel interference (HOMi) 
effect, two photons in otherwise identical modes impinge 
on two different input ports of a 50/50 beam splitter, 
and, thanks to bosonic interference, always emerge to- 
gether in one of the two output ports. More precisely, 
the output state can be expressed in Fock states as 

I*)ab = (|0>a|2>b-|2>aI0>b)/^- (1) 

Here A and B denote the two output modes, with iden- 
tical polarizations, frequencies, and transverse spatial 
quantum numbers, and differing only in their propa- 
gation directions. Great progress has been made re- 
cently in building waveguide circuits on chips, with which 
high-visibility interference fringes involving multi-photon 
states with high purity such as l^*) can be observed [S]. 

The aspect of the output state l^')^^ that interests us 
here is that it, provided the modes A and B are spatially 
separated, is entangled. For instance, the pure state \'^) 
can be shown to violate Bell-type inequalities [6,. What 
concerns us in particular, is how one could verify the 
entanglement of noisy versions of the ideal state, con- 
taining, e.g., phase noise and contaminations with states 
with different numbers of photons (no photons at all, one 
photon in total, or more than two photons in total). As it 
turns out, standard measurements and operations used 
in, e.g., to characterize and manipulate few-photon 
states are indeed sufficient for entanglement verification, 
provided (but this is a far from trivial proviso) all photo 
detectors detect photons only in particular modes. That 
is, if we assume our detectors are sensitive only to one 
particular polarization, spectral profile, and transverse 
spatial mode, then the method we present here will un- 
ambiguously detect entanglement even if the actual in- 
put state (with arbitrary numbers of photons in it) has 
a multi-mode character. Moreover, in this case we can 
construct lower bounds on the amount of entanglement 



as well. The reason is, that such a detection scheme is 
equivalent to a protocol where a filtering operation is ap- 
plied to the input state that keeps only photons in the 
desired modes. Since this operation is local, the amount 
of entanglement of the resulting filtered state cannot be 
larger, on average, than the entanglement present in the 
input state. 

On the other hand, if we drop the assumption about 
the single-mode character of our detection devices, then 
the problem of verifying entanglement of a delocalized 
two-photon state becomes much more involved, also when 
compared to the similar problem of verifying entangle- 
ment of a delocalized single photon [71 18] . We will give 
the essential reason for this difference and present solu- 
tions for the multi-mode multi-photon entanglement ver- 
ification problem that will work if the state under inves- 
tigation is sufficiently close to a single-mode entangled 
state. 

It may be interesting to compare our entanglement ver- 
ification scheme to a scheme proposed in Refs. [9l [TO] . 
which likewise uses the HOM interference effect (but in 
its fermionic version) to detect entanglement. The lat- 
ter scheme detects entanglement between electrons, and 
assumes the number of electrons in each input port of a 
50/50 beamsplitter is fixed and known, whereas we do 
not assume a fixed photon number. Indeed, such an as- 
sumption is perfectly fine for first-quantized electrons, 
but not for second-quantized photons. Moreover, we use 
the inverse HOM effect to detect entanglement in a state: 
ideally, we have either two photons or no photons in each 
input mode, whereas Refs. [TU] consider, in the ideal 
case, one electron in each input mode, and then use the 
proper HOM effect for entanglement detection. 

Finally, we recall that the (proper) HOM effect has 
been used to detect entanglement between two input pho- 
tons (see, e.g., Ref. [IT] and references therein). It's still 
true that the assumption that there is exactly one photon 
in each input port is not warranted in general, but, for 
entanglement verification, it is an allowed filtering oper- 
ation, as it is local. In contrast, filtering on having two 
photons in total in the two input ports (which operation 
we would like to perform for our case) would be nonlo- 
cal. Also note that in our case, the output of the inverse 
HOMi experiment would ideally be a product state of 
two photons. 
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II. ENTANGLEMENT VERIFICATION FOR 
SINGLE-MODE STATES 

A. Defining "single mode" 

Let us first consider so-called single mode states, by 
which we mean states where any photons present are in 
the same transverse spatial, spectral, and polarization 
modes, with the understanding that they can differ in 
their direction of propagation (there are two such modes 
in our case, spatially separated, which we call modes A 
and B). Since experiments typically must be repeated 
in time, we do allow the spectral mode functions (/>(w) 
to differ by a phase factor exp^iujT) with T a known 
delay time, without the photons losing their single-mode 
character. 

We could, in principle, perform tomography on the full 
state to determine its density matrix and from this calcu- 
late a measure of entanglement, e.g., the concurrence or 
negativity of the state, and thus determine whether the 
state is entangled. However, since we shouldn't assume 
anything about the Hilbert space that the state lives in 
(since we want be able to verify the entanglement on 
noisy versions of our ideal state), we would have an infi- 
nite number of matrix elements to determine. Even if we 
were to make restrictive assumptions about the Hilbert 
space of the state, it would still require numerous mea- 
surements to fully determine the state. For example, if 
we assumed that the state did not contain more than 
two photons, this would still leave a 6x6 density matrix 
to determine. If we are not interested in fully character- 
izing the state, but merely in verifying its entanglement 
we do not need to do so much work. Instead of trying to 
exactly calculate a measure of entanglement of the state, 
we can instead calculate a lower bound which will allow 
verification of entanglement of the state with far fewer 
measurements . 



B. Local filtering 

Let the state whose entanglement we are trying to ver- 
ify be called p. A bound on the entanglement can be 
found in the following way. Suppose we were to apply the 
following local filtering operations: we ask about each of 
the two spatially separated modes A and B two questions 

Filter "l" : Is there exactly 1 photon in the mode? 
Filter "2" : Are there more than 2 photons in the mode? 

We consider this filtering a success if the answer is "no" 
to both questions [cf. Eq. ([I])]. The probability then 
of successful filtering is P = Po,o + Po,2 + ^2,o + ^2,2j 
where Pij is the probability to find i photons in mode A 
and j photons in mode B in the unfiltered state p. This 
filtering collapses our state to one living in the smaller 
Hilbert space spanned by |0)^|0)b, \0}a\2)b, \2)a\0)b, 
and |2)yi|2)B. At this point we have a state represented 



by a density matrix with up to 16 nonzero elements. To 
simplify calculations we can further bound the state's 
entanglement by assuming we apply another local opera- 
tion, which in addition requires classical communication: 

Local operation + CC : "phaseshift" : 

apply the same random phase shift to both modes 

thus destroying any coherence between states with dif- 
ferent numbers of photons and reducing the number of 
nonzero matrix elements to at most 6. 



C. Entanglement criterion 



The end result of filtering is of the (normalized) form 
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Since concurrence is an entanglement monotone and p is 
the result of only local operations and classical commu- 
nication applied to p, the concurrence of p bounds the 
concurrence of p: PC{p) < C{p). The concurrence of p 
is 

PC{p) = max[ 0, 2\d\ - 2/fV^ ] (3) 
which is greater than zero when 



(4) 



Thus p is provably entangled if ineq. Q holds true, and 
so too is p. Similarly, since negativity is also an entangle- 
ment monotone, the negativity of p bounds the negativity 
of p in the same way: PN{p) < J^{p)- But calculating 
the negativity of p results in exactly the same bound as 
found by calculating the concurrence: the state is prov- 
ably entangled if Po-P2,2 < MP- 

Now we must find a way to bound Since d — 

P(02|p|20) = (02|p|20) we don't need to physically per- 
form any of the filtering mentioned above, as we can de- 
termine the needed information, d, from the unfiltered 
state p. To do this, consider placing the two modes of p 
on the two input ports of a lossless 50/50 beamsplitter. 
We will label the input modes A and B, and the output 
modes C and D. The transformation between input mode 
creation operators and output creation operators can be 
written as follows (after adding, for convenience, a 7r/2 
phase shift to mode D to compensate for the tt/2 phase 
shift upon refiection) 



+ fit 



and 



(5) 



which allows us to calculate photo-detection probabilities 
Qij for the output modes, where Qt^j is the probability 
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Take, for instance, the state 



Entanglement 




0.01 0.02 0.03 0.04 



FIG. 1: Scatter plot of the right-hand side vs the left-hand 
side of our entanglement criterion Q. Red dots lie on the 
boundary of entanglement vs separable, and correspond to 
pure separable states of the form (|0)a + o|2)>i) <8> (|0)b + 
6|2)_b) where a and h are real. Blue triangles corresponds to 
mixtures of two randomly generated separable states of the 
form (|0)a + ai|l)A + a2|2)A) ® (|0)s 6i|l)s 62|2)s) (with 
complex coefficients). 



to find i photons in mode C and j photons in mode D. 
It can be shown that 



(6) 



which gives 



So when 



-P2.0 + -Fb,: 



(7) 



(8) 



the state can be said to be provably entangled. Figure 1 
plots both sides of our inequality ([8| for many randomly 
picked separable states, to show how this criterion indeed 
verifies entanglement. Moreover, the figure caption iden- 
tifies the states lying on the borderline between separable 
and verifiably entangled. 



D. An additional phase shift 



pi := ^ |00) (00| + i(|20)+z |02))((20|-z (02|)+^ |22) (22| . 

(9) 



For this state \d\'- 



and Po-P2,2 = 



36 



so by ineq. 



(|4|) the state is in fact entangled. However 3?(d)^ — 0, 
so ineq. (Sj) will not detect the entanglement. But if we 
apply a phase shift of exp(i^) to one of the modes then d 
will become purely real (and so (5i,i will be maximized), 
and ineq. (|8]) will detect the entanglement. As Figure 2 
(top) shows, for this state with a phase exp(i(/()) applied 
to the first mode, entanglement will be detected when 

TT. A similar. 



is between \'k and Itt or between Itt and 

6 6 6 

but more noisy state. 



P2 := 3 |00) (00K-(|20)-fz |02))((20|-i (02|)+- |22) (22| , 

(10) 

will have a smaller range of detectable entanglement, 
specifically when is between .397r and .GItt or between 
1.397r and I.GItt (see Figure 2, bottom part). 





Our condition ([8]) will not detect entanglement in an 
input state, even when it is in fact present, when d is 
largely or purely imaginary. But if one were to apply a 
phase shift to one of the modes before placing the state 
on the beam splitter and vary that phase until Q\^\ was 
maximized (the same local operation with classical com- 
munication as performed in [5 ), this would maximize 
3fi(d)^, thus making ineq. ([4]) equivalent to ([8|. In other 
words, such states then can be detected by our criterion. 



FIG. 2: (Qi,i - |(P2,o + Pq,i)) , that is, the right-hand side 
of inequality ([sjl, for the state p\ (top), defined in (|9|, and 
the more noisy p2 (bottom), defined in (10 1, as a function of 
a phase shift exp(i0) applied to the first mode. The shaded 
region represents for which values of entanglement will be 
detected by ineq. (|8| [both states are entangled for any value 
of (Al. 
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E. Asymmetric beamsplitters 

To bound d we placed our state on a 50/50 beamsplit- 
ter, but it is easy to generalize our analysis to beam spit- 
ters which are not equally balanced. Suppose our beam 
splitter has a (real) reflection coefficient r and a (real) 
transmission coeSicient t = y/l — such that the input 
creation operators transform as 



-> rc^ + td) 



and 



tc^ - rd'^ 



(11) 



Following the same analysis as before we find that if 

P^P2a<y 2 ) 

the state is provably entangled. 



III. CONSIDERATIONS CONCERNING 
MULTI-MODE MULTI-PHOTON STATES 

We made the assumption at the beginning of our anal- 
ysis that any photons present are in the same transverse 
spatial, spectral, and polarization mode. However if our 
detectors only detect a certain single mode we can drop 
the assumption that the photons being in the same mode 
as this is equivalent to a local filtering. That is, using sin- 
gle mode detectors is equivalent to an additional filtering 
performed on each of the spatially separated modes, fil- 
tering out all photons not in the single mode of interest 
before detection takes place. What if we drop the as- 
sumption of single-mode detectors? 

Suppose we have an input state in which the photons 
present are not all in the same transverse spatial, spec- 
tral, and polarization mode. The entanglement verifi- 
cation scheme described above did assume that the two 
photons in the filtered state (after the local filtering op- 
erations "1" and "2") are in the same mode, because of 
the explicit assumption that there is interference (of the 
"inverse HOM" type) taking place on a beam splitter. 
But this assumption does affect how we interpret the re- 
sults of the measurements: in particular, the quantity 
Qii (which we would like to be large) could be danger- 
ously contaminated with contributions from those input 
states that lead to larger values of Qn for photons in 
different modes than for photons in the same modes. For 
example, if we start with an output state with one photon 
in each output port, but of different colors, then applying 
the inverse beam-splitter transformation yields an input 
state that has this undesired property. The question is 
to what extent we can avoid or correct for the presence 
of such input states. 



of state, i.e., one that maximizes the right hand side of 
Eq. ([s]) without HOM entanglement, such as the state 
mentioned above 



(|10).e<i+|01).ed)®(|10)i 



|01)Mue)/2 (13) 



While this state has twice as much entanglement as the 
HOM state, it is not the type of entanglement we are 
interested in trying to detect here. A state such as this 
with a probability P2 of detecting two photons of differ- 
ent color will contribute at most ?>P°i/2 to the quantity 
being squared on the rhs of Eq. Q , so we will compen- 
sate for this possible contribution by subtracting 
For states close to the ideal state the contamination of 
different colored photons will be small and thus the cor- 
rection will be small. We can also construct a bound that 
does not rely on measuring the probability of detecting 
two photons of different colors, since it is always less than 
or equal to the probability of detecting two photons of 
any co\oi{P° ^ < Pi^i.) Using this, our (conservative) 
condition for entanglement becomes 



^0^2,2 < (niax[ Qi,i - Pi,i - P2/2 , 1)^ 



(14) 



B. Nonexistence of local filters for sameness of 
modes 

It would be nice if we could find a local filtering opera- 
tion that checks whether two input photons propagating 
in one direction are in the same mode with respect to 
the other quantum numbers or not. There is certainly 
no von Neumann measurement that achieves that goal, 
as the target states are not all orthogonal. But, surpris- 
ingly, we cannot even construct a POVM that does the 
trick: the reason is that even if we start with a state that 
contains two photons in orthogonal modes, say described 
by creation operators a\ and 03, then we can view the 
same state as a superposition of two states, each with 
the two photons in identical modes, as described by the 
creation operators aj_ — (a| ±a2)/\/2. This results from 
the identity 



t t («l)' - («-) 



-.t \2 



(15) 



This is then the essential difference between single- 
photon states and multi-photon states, which makes 
entanglement verification much harder for two-photon 
states than for single-photon states! Moreover, this also 
illustrates a difference between bosons and fermions: in 
the case of two fermions there is an antisymmetric sub- 
space, and, e.g., we can certainly perform a measurement 
that checks whether two spin- 1/2 systems have different 
spins (singlet state!) or not. 



A. A corrected entanglement criterion 



C. An alternative local operation 



One way of correcting for these unwanted states is to 
subtract the contribution from the worst possible kind 



All is not quite lost, as we can still apply other sorts of 
local operations that are useful for the analysis of entan- 
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glement of the input state. In particular, suppose that 
our input state is some coherent superposition of, e.g., 
the desired state |0)^ |2)g — |2)^ \0)g and an unwanted 
state |1)a2 (with photons in different modes). 
There is a local operation that transforms this superpo- 
sition into an incoherent mixture of these two states: for 
each pair of orthogonal modes Ak and Bk (picked from 
some fixed basis: that's the essential difference from the 
no-go statement from the preceding subsection) apply 
a random fc-dependent phase shift, and then forget the 
precise phase shifts applied. This operation will only pre- 
serve the coherence of superpositions of photons in the 
same spectral, polarization and transverse modes in A 
and B. That is, by a local operation we can transform 
the input state into a state of the form 



PsP' 



a -p. 



(16) 



where the first term denotes states that do display (in- 
verse) HOM interference, and the second term states that 
do not; Pg is the probability of observing HOM interfer- 
ence, given p. The point is that we have now separated 
the input state in two parts, the first part of which is 
the state for which our method demonstrates entangle- 
ment (see below for further elaborations of this point). 
The second term has no entanglement, since any super- 
positions in that term have been destroyed. Its presence 
could imply the state p is not entangled, even if p'^ is, 
namely if 1 — Ps is too large. We will not solve the (hard) 
general problem of identifying for what values of Ps and 
for what states p*, entanglement of the latter still implies 
entanglement of p. 

Let us return to the statement that is entangled, if 
our verification method succeeded. We still have to dis- 
cuss the fact that our method assumed that both pho- 
tons are in one particular mode, whereas for photons in 
we only know they are in the same mode, but not in 
which one. This does have consequences for the amount 
of entanglement (see [12] for extensive discussions of this 
issue), but not for the bare fact that the state is entan- 
gled. We can demonstrate this by showing that the state 
p'* can be distilled (the following protocol is far from op- 
timal, and one can easily improve its efficiency; here its 
point is only an existence proof): just take two copies of 
p"; first determine a particular mode such that the pro- 
jection of p* onto that mode is entangled; then perform 



on each of the A and B modes a joint measurement that 
counts how many photons in that particular mode there 
are in total in the two copies. If the answer is "2" for 
both A and B, we have an entangled state in that one 
particular mode. In this highly inefficient protocol the 
average amount of entanglement decreases (unless only a 
single mode is occupied), but it stays nonzero. Hence p^ 
must be entangled. 

IV. SUMMARY 

We demonstrated how the inverse HOMi effect can be 
used to verify the mode entanglement present in a state 
of the form (|0) |2) — |2) |0))/\/2, and noisy versions of 
it. If the photons in the state are all "single- mode" , 
that is, all have the same polarization, the same trans- 
verse mode profile and the same spectral profile, then 
our method easily bounds the amount of entanglement 
from below. That directly gives a criterion, inequality 
([S]), which, when satisfied for a given single-mode state, 
is sufficient to prove entanglement. We analyzed how 
the applicability of the criterion can be improved simply 
by applying an additional phase shift to one of the two 
modes. The operations needed to verify entanglement 
can be implemented with linear optics, and are just those 
demonstrated in the experiment of ;5]. 

We discussed how the problem of verifying entangle- 
ment in the delocalized two-photon state with the inverse 
HOMi effect becomes more "interesting" (a euphemism 
for "complicated" ) without this single-mode assumption 
[more precisely, when both the input state and one's 
photo detectors are multi-mode], and why a delocalized 
single-photon state does not suffer from these complica- 
tions. On the other hand, the interpretation of violating 
a Bell inequality with unbalanced homodyne measure- 
ments [6; is immune to the single-mode or multi-mode 
character of the input state, at the small cost of requir- 
ing phase-locked local oscillators, thus showing an advan- 
tage of Bell inequalities in the context of entanglement 
verification. 

We gave a simple solution to the full problem of in- 
verse HOMi multi-mode multi-photon mode entangle- 
ment, based on bounding the deviation of the actual 
state from a single-mode state, which works well when 
that deviation is sufficiently small. It yields a (more con- 
servative) entanglement criterion ( 14 ) . 
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